Spin fluctuations in cuprates as the key to high T c 
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Abstract. Spin fluctuations represent the lowest established energy scale in cuprates and are crucial 
for the understanding of anomalous normal state properties and superconductivity in these materials 
(U. The memory-function approach to the spin response in the t-J model is described. Combined 
with numerical results for small systems it is able to explain the anomalous scaling at low doping and 
the crossover to the Fermi-liquid-like behavior in overdoped systems. Within the superconducting 
phase the theory reproduces the resonant peak and its peculiar double dispersion. Such spin fluctu- 
ations are then used as the input for the theory of superconductivity within the t-J model, where we 
show that an important role is played also by the next-nearest-neighbour hopping parameter t 1 . 
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INTRODUCTION 

The phase diagram of cuprates still represents one of the major challenges in solid state 
physics, both for theoreticians and experimentalists. Besides superconductivity (SC) and 
antiferromagnetic (AFM) ordering, several regimes with distinct electronic properties 
have been identified within the normal metallic phase. In this contribution, we are con- 
cerned with the spin dynamical response in cuprates which has been intensively studied 
using the inelastic neutron scattering (INS) Jlllil] and NMR relaxation experiments 0]. 

There is an abundant evidence that in underdoped cuprates magnetic properties are 
not following the usual Fermi-liquid (FL) scenario within the metallic state above the 
SC transition T > T c . It should be reminded that within a normal FL the dynamical 
spin susceptibility %^(cd) is essentially T- independent at low T < Tfl- From the 
point of spin response one can use the latter criterion as the working definition of the 
Fermi-liquid temperature Tfl. Note that in usual metals Tpi ^> 1000 K, in striking 
contrast to underdoped cuprates with Tfl < T c JH 0]- Evidence for the non-Fermi- 
liquid (NFL) behavior are INS results for the q-integrated spin susceptibility which 
exhibit in a broad CO, T range an anomalous, but universal x'K®) 06 /( ft) /^)' established 
in underdoped La 2 xSr x Cu0 4 (LSCO) S H, as well as in YBaCu 3 6+x (YBCO) 
Similar conclusions arise from evident T-dependence of 63 Cu NMR spin-lattice 
relaxation rate \/T\T and of the spin- spin relaxation rate I/T2G M\- 

Cuprates at optimum doping, and even more in the overdoped regime, approach closer 
the usual FL description. INS reveals only weak spin response at low energies co in 
the normal state (NS) at T > T c , characteristic for metals with a broader band. Also, 
NMR relaxation rates \/T\T and I/T2G are weakly T-dependent, again consistent with 
the FL scenario. Analogous message arises from the analysis of cuprates doped with 
nonmagnetic Li and Zn [7], where the impurity-induced spin susceptibility varies as 



\/(T + Tk), whereby the characteristic (Kondo-type) temperature is Tk ~ in the 
underdoped regime and increasing fast with doping in the overdoped regime. 

While in the NS the dynamical spin response %'q((Q) is in general compatible with 
an overdamped collective mode, the prominent feature appearing in the SC phase is the 
magnetic resonant mode. First observed in optimally doped YBCO BSD, it has been in 
the last decade the subject of numerous INS experiments yD. The resonant response, 
in spite of evident differences between YBCO and LSCO systems, as well as changes 
with doping, reveals some surprisingly universal characteristics. The peak intensity is 
highest at the commensurate wavevector Q = (k,k), while its frequency co r increases 
with doping up to the optimum doping. In addition, one component of the resonant mode 
disperses downwards [9], while another branch apparently emerging from the same peak 
shows an upward dispersion 

It seems a plausible (although surprisingly not generally accepted) conclusion, that the 
understanding of spin dynamics is the key for the proper description of the anomalous 
properties of cuprates, and to the mechanism of high T c in particular. The main argument 
remains, that up to now the spin fluctuations represent the lowest (experimentally well 
established) energy scale in cuprates, both in the NS as well as in the SC state. Namely, 
the peak in the spin response (as measured by INS) in the normal state appears in 
underdoped cuprates at co p ~ T c , moving even lower with decreasing T. Moreover, the 
collective mode (O r in the SC phase lies even below the SC gap co r < 2Aq. On the other 
hand, low spin-fluctuation energy scale also sets a clear limit to the FL behavior since 
FL can become normal only for T, CO which are below this scale. 

A comprehensive theoretical description of spin fluctuations in cuprates and their im- 
plications on other properties, in particular their role in the mechanism of SC, is still 
lacking. A FL behavior in the overdoped regime far from a metal-insulator transition 
seems plausible, nevertheless a solid theoretical approach is missing even in this regime. 
A crossover from a strange metal to a coherent metal phase has been predicted within 
some theoretical approach. Quite fashionable and frequently invoked interpretation is 
given in terms of the quantum critical point (QCP) at optimum doping c* h (masked, how- 
ever, at low T by the SC phase), dividing the FL phase at c% > c\ from a (singular) 
non-Fermi-liquid (NFL) metal at c/, < c\. Such a scenario is established, e.g. theoreti- 
cally in spin systems [13] and experimentally in some heavy-fermion compounds, but 
remains controversial in cuprates. The obvious argument against the QCP scenario is 
the absence of a critical length scale, e.g., AFM correlation length £; (T — > 0) — > °° as 
well as the absense of the phase with the AFM long-range order for c/, < c* h (sometimes 
put in connection with the pseudogap scale T*). In the underdoped cuprates the spin 
fluctuation seem to follow quite well the phenomenological scenario of the marginal FL 
ifLUl . which got so far only partially a solid theoretical foundation. The main distinction 
to the QCP scenario is the absense of a critical length scale, which is in agreement with 
low-energy INS revealing at low T the saturation of the inverse AFM correlation length 
K = at least in YBCO £} and in LSCO at low doping flQ]. 

With respect to the most challenging problem, the mechanism of SC in cupratres, 
the role of strong correlations and the antiferromagnetic (AFM) state of the reference 
insulating undoped compound has been recognized very early ffl^l . Still, up to date 
there is no general consensus whether ingredients as embodied within the prototype 
single-band models of strongly correlated electrons are sufficient to explain the onset 



of high T c . Even within the frequently invoked t-J model, being the subject of this 
paper, proposed mechanism of SC and the methods for the evaluation of correponding 
T c differ with respect to the fact whether the attractive interaction is mainly local 
and instantaneous lfl6ll or the retardation effects are important lfl7ll . Recognizing the 
very low spin-fluctuation scale, we will advocate in the following the latter for spin- 
fluctuation scenario, emerging in contrast to previous approaches directly from the 
strongly correlated t-J model. 

In the following we present some of our recent theoretical results on spin fluctuations 
in cuprates and their relation to SC. The analysis within the t-J model, as relevant for 
cuprates, is mostly based on the general memory-function approach and the equations- 
of-motion (EQM) method, The latter has been first applied to the t-J model to explain 
anomalous (MFL-type) properties of NS spectral function [18] and then extended to low- 
doping regime lll9ll and SC lE(ill . Spin dynamical response ^q(co) has been considered 
within analogous treatment to yield the overdamped mode in the NS and resonant 
peak dispersion in the SC state lElll . the anomalous oj/T scaling in the underdoped 
regime lE^lhthe influence of nonmagnetic impurities ll23ll, th e NFL-FL crossover in spin 
dynamics II24I1 . and double dispersion of resonant peak »2jj\ . The extracted knowledge 
on spin fluctuations is used as an input the theory of SC 12611 . 



NFL - FL CROSSOVER IN THE NORMAL STATE 

To be specific, we consider in the following the spin dynamics within the framework 
of the extended t-J model, which has been shown to represent surprisingly well several 
electronic properties of cuprates, both qualitatively and quantitatively 12711 . 

H = - £ Ujc] s Cig+Jj^{Si -Sj- -mrij), (1) 

{ij) 

including in general both the NN hopping ttj = t and the NNN hopping t\j = t', and 
involving the projected fermionic operators, Si s = (1 — n,i- s )ci s . 

We will first argue 1E2I1 that the anomalous oj/T scaling and the related NFL behavior 
of the magnetic response can be understood as a consequence of few simple ingredients 
which appear to be valid for doped AFM in the normal state: a) the collective mode is 
strongly overdamped, whereby the damping is nearly CO- and T- independent at low GO, 
and b) there is no long-range spin order at low T, so that static spin correlations saturate 
with a finite £ . 

Within the memory function approach the dynamical spin susceptibility Xq(®) = 
— ((Sq;5q))a> can be generally expressed |2jJ,|22|] in the form 

Zq(w) = (& + ojM^co) _ col > ^ 

suitable for the analysis of the magnetic response, as present in undoped and moderately 
doped AFM [21]. (0 q represents the frequency of the collective mode provided that the 
mode damping is small, i.e., A q ~ Mq(ft)q) < co q . In the opposite case, i.e. A q > co q the 



mode is overdamped. Still, the advantage of the form © is that it can fulfill basic sum 
rules even for an approximate M q . Thermodynamic quantitites entering Eq. © can be 
expressed as 

nq = -i([Si q ,^)]), < = V*2> (3) 

where Xq = Xqi® = 0) is the static susceptibility. 

7] q is the spin stiffness and can be expressend in terms of the static correlation 
functions, in particular within t-J model t]q = —(H)/N. Xq ( or w q) remains to be 
determined, even for known M q (co). It is quite a sensitive quantity, hence it safer to 
fix it by the sum rule (fluctuation - dissipation relation) 

1 f°° CO 

-j Q dco cth— 4(co) = (5i q 5 q )=C q , (4) 

given in terms of equal time spin correlations, which are expected to be much more 
robust. Moreover C q are bound by the constraint (1/Af) £ q C q = (1 — Q,)/4, where c/, is 
an effective hole doping. 

Let us now state two basic assumptions: a) static correlations are taken to follow a 
Lorentzian form, i.e. C q = C/(K 2 + q 2 ) where q = q — Q. Kis assumed to be a noncritical 
quantity, which on approaching low T saturates at a finite value. As already noted, this 
is consistent with the neutron scattering data for weakly doped LSCO o] and YBCO 
It is also consistent with numerical results for the t-J model at finite doping, b) 
The damping is also assumed to be a constant, M q (co) ~ A, i.e., (roughly) independent 
of CO , q and T, or at least not critically dependent on these variables. The support for 
this assumption comes from our numerical results on small systems, using the finite- 
T Lanczos method (FTLM) ll27ll . for the t-J model on small lattices with up to 20 
sites. Calculating Z q (co) and extracting then M q (co) with the help of Eq.©, one can 
conclude |22t| that in spite of widely different Xqi®) the damping function M q (co) is 
nearly constant in a broad range of co < t and almost independent of q. Moreover, for 
doped systems with Ch > data are consistent with a finite (and quite large) extrapolated 
value Aq(T — > 0) . In the normal state this leads to a overdamped collective mode vicinity 
of q = Q, i.e., COq < A, as generally observed in INS experiments ASS]* 

and r q < COq. At low T — > 0, Eq.© leads now to a nontrivial restriction for co q and r q . 
The relevant quantity is the peak frequency CO p = Tq(T — » 0), which determines the 
characteristic T = spin-fluctuation scale as well as Tfl- 
The crucial parameter appears to be 

C =CTzA/(2r)K 2 ), G) p ~Ae~ 2? , (6) 

which exponentially renormalizes co p . Since C ~ 0(1) and r\ ~ 0.6 t at low doping, £ is 
effectively governed by the ratio A/k 2 . It is easy to imagine the situation that £ 3> 1 in 
the underdoped cuprates, leading to very low C0„ A and even co„ < T c . On the other 



hand, in the overdoped case A/zc 2 ~ 1 and (O p becomes large as in usual FL systems. 
Due to exponential dependence in Eq.© it is also plausible that the crossover from the 
NFL regime with extremely small (O p and FL behavior is quite abrupt 12411 . resembling 
the QCP scenario. 

In order to extract the characteristic energy scale (Ofl of spin fluctuations directly 
from numerical FTLM results Il24ll . we use an alternative definition, 



<Qfl(T)= S q /Xq(T) 

with the corresponding T = limit (Ofl(0). Note that (Ofl(0) 
frequency moment of the shape function %q(co) / CO, 



(7) 

(to) is just the first 



co FL (0) = ((o) 



Jo 



XQ(co)dco. 



(8) 



On the other hand, one can extract (Ofl also from experiments, in particular from NMR 
1/?2g relaxation data [4], which give rather straightforward information on %q(T). 

In Fig. 1 we show FTLM results for <»fl(qO at O.lt < T < J. Besides we also present 
values extrapolated to T — > 0. Note that in the considered T window Sq(T) is essentially 
T-independent, following well the linear variation 1/Sq = Kc%. In contrast, the FL scale 
(Ofl reveals a nonuniform variation with doping. Again, for q, > c* h , (Ofl is already 
rather T-independent for T < J. On the other hand, in the regime c/, < c* h we find a 
strong T-dependence of (Ofl even at lowest reliable T, where (Ofl ~ T + 0)fl(0) • We can 
summarize results in Fig. 1 as follows: a) in the overdoped regime (Ofl(0) ~ a(c/, — c/,o) 
with CM ~ 0.12 and a large slope a ~ 3.5? ~ 1.4 eV, b) in the underdoped regime our 
results indicate on a smooth crossover to very small (Ofl(0) <C J- 
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FIGURE 1. FL scale ©fx/f vs. c/,, obtained for the t-J model using the FTLM for T > and their 
T — > extrapolated values. The inset shows T — results for 1 /Sq vs. c/ 7 . Dashed lines are guide to the 
eye only. 



Let us estimate Xq(T) and consequently (Ofl directly from experiments on cuprates. 
Within the normal state the results for the NMR spin-spin relaxation time Tig, obtained 



from the 63 Cu spin-echo decay, can be related to static Xq 112811 . Assuming that Xq is 
peaked at commensurate q = Q and can be described by a Lorentzian form with a width 
K <tLK, one gets a simplified relation 



J-~0.083kF(Q)^ q . 

t2G 



(9) 



I/T2G relaxation rates have been measured and summarized in Ref. 0], i.e., from 
underdoped to optimally doped YBCO with 0.63 < x < 1, underdoped YE^CuztOs, 
nearly optimum doped Tl2Ba2Ca2Cu30io (Tl-2223) and the overdoped Ti2Ba2Cu0 6+ g 
(Tl-2201), whereby the normalization with corresponding F(Q) has been already taken 
into account. Relevant K is the one appropriate for low- ft) spin dynamics and measured 
directly by INS. For YBCO data are taken from Ref.[29], which allows us to evaluate 
Xq(T) from Eq.©. Sq is so far not experimentally accessible, so we assume here the t-J 
model results to finally extract corresponding (Qfl{T) as presented in Fig. 2 for various 
cuprates. Derived O)fl(O) are well in agreement with model result in Fig. 1, in particular 
regarding the large slope in the overdoped regime and a clear change of scale between 
the underdoped and overdoped cuprates. 



200 



150 



> 

— 100 



50 




o YBCO c 



100 200 300 400 
T[K] 



500 



FIGURE 2. (£>fl vs. T, for various cuprates. The inset shows the extrapolated scales Ofx(0) and © 
(defined by I/Xq 05 T + ©) vs. doping c/ 7 . 



DISPERSION OF THE RESONANT MODE 

Using the method of equations of motion within the t-J model it has been shown that 
the collective spin fluctuations decay into electron-hole excitations 112 lL 12511 . This leads 
to the lowest-order mode-coupling approximation for the damping in the NS, 

Aq(») = ^^/^1/(«0-/(« + «0]E^k(^k+q(» + W / ), (10) 



where Wkq is the effective spin-fermion coupling [21] and A k (co) is the single-particle 
spectral function. Provided the existence of 'hot spots' where the FS crosses the AFM 
zone boundary (being the case for cuprates at low to intermediate doping) we assume 
that in the NS low- ft) quasiparticles (QP) with dispersion £ k and weight Zk determine 
the spectral function Ak( ft)) = Z k 8((0 — £k). This results in a rather constant A q ( ft)). The 
form of Eq. (flOb is anyhow quite generic for the damping of the collective magnetic 
mode in a metallic system, since the lowest-energy decay processes naturally involve 
the electron-hole excitations close to the FS. Similar expressions appear also in theories 
based on the RPA approach ll30l l3lll. Within the SC phase, Eq. (flOb has to be generalized 
to include the anomalous spectral functions [30] leading to 12 ll 12511 

Aq(ft>) ~ 2^vE # kq( M kVk+q - V k Wk+q) 2 [/(^k) ~ f(E k ~ ft))] 8 (CO - E k - E k+q )] , 

(ID 

where w^ q = w^ q ZkZk+ q /7]q, while Wk,Vk are the usual BCS coherence amplitudes and 



+ A£. For the SC gap we assume the d x 2_ y i form, A q = Ao(cos^ x — cosq y )/2. 
Thus we end up with few adjustable parameters at chosen c/ 7 : K in the Lorentzian form 



for C q , the effective coupling w and the maximum SC gap Aq 112 lL 12511. 





FIGURE 3. Xq (®) at intermediate doping c/, = 0.15 for momenta: a) along the x direction q = ^(1,0), 
and b) along the zone diagonal q = q(l, 1). 

At intermediate doping the collective mode is heavily overdamped in the NS. The 



indication for the latter is low intensity of the INS in the relevant low-energy window. 
For the presented case 12511 we fix the 'optimum' doping at = 0. 15 and k ~ 1 .25. The 
SC gap is roughly known from experiments Ao = 40 meV. The remaining input is Tq 
within the NS. For the appearance of the upper resonant branch it is crucial that Tq is 
not too large, as seems to be inherent within the RPA II30L I3 ill. 

In Fig. 3 we display Y.!1{(0) for momenta both along the two q directions. Following 
observations can be made 12511 : a) results reveal two branches emerging from the same 
coherent resonant mode at co, ~ 41 meV. Intensity plots of both branches within the 
q plane are square-like around AFM Q, however with quite pronounced anisotropy. b) 
For the downward branch the intensities are strongest along the (1,0) direction, c) The 
development is more sensitive for co > C0 r , still the situation with the upward branch is 
just opposite to the downward one. The dispersion is stronger along the (1,0) direction, 
d) Above the damping threshold CO > 2Aq the upward branch merges into an incoherent 
response broad both in q as well as in CO. The incoherent part still exhausts most of the 
intensity sum rule, Eq. ©, even for q = Q. 



SPIN-FLUCTUATION MECHANISM OF SUPERCONDUCTIVITY 

The projection in fermionic operators in the t-J model, Eq.®, leads to a nontrivial EQM, 
which can be in the k-basis written as 

e° 1 1 

[Cfo, H] = [(1 +C h ) -± - (1 - C h )J]£te + -j= V m kq [^C k _ q ^ + Sjcfc q. j - -HqCk-qJ , 

1 Viy q l 

(12) 

where m kq = 27 y q + is the effective spin-fermion coupling, while = —Aty^ — 
4t'Y k is the bare band dispersion on a square lattice, and y q = (cosk x + cosk y )/2, Y q — 
cosk x cosk y . To keep similarity with the spin-fermion phenomenology lfl7ll we use the 
symmetrized coupling lElll 

m kq = 277q + (e0_ q + e0)/2. (13) 

EQM, Eq. (fT2"l) . are used to derive the approximation for the Green's function (GF) 
matrix G ks (co) = ((Vbl^L)) o for the spinor ^ = (c k)5 ,cl k A We follow the 
method, as applied to the normal state (NS) GF by present authors HHE^], and gener- 
alized to the SC pairing in Ref.JHEl. In general, we can represent the GF matrix in 



the form 

G k ,(w) 1 = -[coT -^ s + ^T 3 -L ks ((o)}, (14) 

LA/ 

where a = Y,i({cis, c] s }+) /N = (1 +c^)/2 is the normalization factor, ji is the chemical 

/V J. 

potential, and the frequency matrix £ kl . = ({[ v I / ks ,//], v I'^ 1 ,} + )/a, which generates a 
renormalized band £ k = ^ = £ — Arjity^ — 4r]2? / 7 k and the mean-field SC gap 

= £L 2 = -^E^-q^-q,-^)- d5) 



To evaluate £^({0) we use the lowest-order mode-coupling approximation, analogous 
to the treatment of the SC in the spin-fermion model 11711 . Taking into account EQM, 
Eq. (fT2"l) . and by decoupling fermionic and bosonic degrees of freedom, one gets 

^11(12),. \ — 3 y-, _ 2 ^H(12)/. \ /• \ /i^x 

E k, ( l( °n) = L m kq G k qj (^)*q( Z£0 « ~ lG) >n) ( 16 ) 

where i(O n = in(2n+ l)//3 and /3 = l/T, whereby we have neglected the charge- 
fluctuation contribution. 

In order to analyze the low-energy behavior in the SC state, we use the QP approx- 
imation for the spectral function matrix where the QP energies are = (e^ + AjL) 1 / 2 , 
while NS parameters, i.e., the QP weight Zk and the QP energy £k, are determined from 
Gks((0 ~ 0), Eq. (fT4ll . By defining normalized F q (iC0i) = %q{}(Oi) / an d rewriting the 
MF gap, Eq. (fT5t . in terms of the spectral function, we can display the gap equation in a 
BCS-like form form, 

^ = ^D 4 ^ qZtqCq.k q]^^th(^), (17) 

where Ck q = Ikq(i6>n ~ 0)//^ pl a Y s the role of the cutoff function with 

and/0 = th(/3£ k /2)/(2£ k ). 

Eq. (fT7t represents the BCS-like expression which we use to evaluate T c . To proceed 
we need the input of two kinds: a) ^q(d)), and b) the NS QP properties Zk,£k- As 
discussed in Sec. II the NS spin dynamics at q ~ Q is generally overdamped in the 
whole doping regime [3]. Hence we use the form, Eq.©. We end up with parameters 
ZqjTq, fc, which are dependent on Ch, but in general as well vary with T. Although one 
can attempt to calculate them as described in Sec. Ill [21], we use here the experimental 
input for cuprates, as discussed in Sec. II [24]. For the NS A^(co) and corresponding 
^k, £k we solve Eq. (TToT) for E k ! = Ek as in Ref. Ill 911 . with the same input for Xq((o). The 
main message remains [26] that soft AFM fluctuations with q ~ Q lead through Eq. (fT6l) 
to a reduction of Zk, which is k-dependent. A pseudogap appears along the AFM zone 
boundary and the FS is effectively truncated in the underdoped regime with Zk F 1 
near the saddle points (n, 0) (in the antinodal part of the FS) I I 9i 12611 . 

Close to half-filling and for %q peaked at q ~ Q both terms in the gap equation, 

Eq. (fTTt . favor the d x i_ y i SC. The mean-field part A^, Eq. (fT5l) . involves only J which 
induces a nonretarded local attraction, playing the major role in the RVB-type theories 
II15L [la]. In contrast, the spin-fluctuation part represents a retarded interaction due 
to the cutoff function Ck q , determined by I\ q . The largest contribution to the SC 
pairing naturally arises from the antinodal part of the Fermi surface. Meanwhile, in the 
same region also Zk is smallest thus reducing the pairing strength, in particular in the 
underdoped regime. 



One can question the relative role of the hopping parameters t,t' and the exchange J 
in the coupling, Eq. ( TT31 ). While our derivation within the t-J model is straightforward, 
an analogous analysis within the Hubbard model using the projections to the lower and 
the upper Hubbard band, respectively, would not yield the J term within the lowest order 
since J t 1 . This stimulates us to investigate in the following also separately the role of 
J term in Eq. (flTl) . both through the MF term, Eq. (fT5t . and the coupling mtq, Eq. (fj~3b . 
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FIGURE 4. T c /t vs. doping q, for t'/t — —0.3, calculated for various versions of Eq. (I17> : a) full result 
(full line), b) with neglected MF term (dashed line), and c) in addition to b) modified m^q without the J 
term (dotted line). 
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FIGURE 5. T c /t vs. —t'/t for fixed 'optimum' doping c/, = 0.17 and different versions of Eq. (I17> . as 
in Fig. 2. 

Results for the NS spectral properties reveal that the coupling to AFM fluctuations 
partly change the shape of the Fermi surface, more pronounced is however the effect on 
the QP weight. is reduced along the AFM zone boundary away from the nodal points. 
In Fig. 4 we present final results for T c , as they follow from the gap equation, Eq. (fTTb . 
It is evident that the spin-fluctuation contribution is dominant over the mean-field term. 



When discussing the role of the J term in the coupling, Eq. (TH51) . we note that in the 
most relevant region, i.e., along the AFM zone boundary m^g = 27 — At' cos k x . Thus, 
for hole doped cuprates, t' < term and J term enhance each other in the coupling, and 
neglecting J in reduces T c , although at the same time relevant is enhanced. 
Finally, in Fig. 5 we present results, as obtained for fixed intermediate doping Cf, = 

0. 17. but different t'/t < 0, as relevant for hole-doped cuprates [32]. As expected the 
dependence on t' is pronounced, since the latter enters directly the coupling m^Q, 
Eq. (fT3l . It is instructive to find an approximate BCS-like formula which simulates our 
results. The latter involves the characteristic cut-off energy Tq, while other relevant 
quantities are the electron density of states JVq and Z m being the minimum Zk on the FS 
(in the antinodal point). Then, we get a reasonable fit to our numerical results with the 
expression, 

T c ~o.5r Q e- 2/( -^ Ve ff\ (19) 

where the effective interaction is given by V e ff = 3Z m (27 — 4t') 2 XQ- 

Probably the most interesting novel result on SC is a pronounced dependence of T c 
on t' which is also consistent with the evidence from different families of cuprates [32]. 
One can give a plausible explanation of this effect. In contrast to NN hopping t, the NNN 
t' represents the hopping within the same AFM sublattice, consequently in a double unit 
cell fermions couple directly to low-frequency AFM paramagnons. Calculated T c are in 
a reasonable range of values in cuprates. We also note that rather modest 'optimum' T c 
value within presented spin-fluctuation scenario emerge due to two competing effects 
in Eqs. (fT7t.(IT9l): large and Xq enhance pairing, while at the same time through a 
reduced Z^ and cutoff Tq they limit T c . At the same time, INS experiments [3] reveal 
that in underdoped cuprates the resonant peak at co ~ (O r takes the dominant part of 
intensity of q ~ Q mode which becomes underdamped possibly even for T > T c . Thus it 
is tempting to relate Tq to (O r and to claim T c ~ a(O r , as indeed observed in cuprates 01 . 
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